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Abstract 

The application of the two-photon process pp ^ pp + for the luminosity mea- 

surements at LHC with the ATLAS detector is considered. The expected accuracy 
of the absolute offline luminosity determination is f H- 2 % for the luminosity range 
of 10^^ 10^^ cm~^s~^. The preliminary cross section estimates done for LHCb 
promise the same level of the luminosity measurement accuracy at L = 2 • 10^^ 

— 2 —1 

cm s . 
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1 Introduction 



The possibility to use the two-photon pair production for luminosity mea- 
surements at hadron colliders was first considered in [1]. It was shown that 
the cross section of the process pp ^ pp + e~^e~ for the forward region can 
be calculated within QED with an accuracy of 1% or better allowing for a 
precise luminosity determination using the number of events observed and the 
calculated cross section: L = N°^'^/a°^'^. 

The method suggested in [1] provides for both the high rate monitoring and the 
absolute luminosity measurements but requires a dedicated apparatus allowing 
for the e+e~-pair detection in the pseudorapidity range of \rj\ > 5 and capable 
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of the pt measurements with the accuracy of a few MeV. This approach has 
never been used. The kinematic domain of the central rapidity pairs is much 
more attractive from the experimental point of view though the corresponding 
cross section is smaller and some additional theoretic uncertainties appear. 

The application of the centrally produced e+e~-pairs with = 0.4 -h 1 GeV 
for luminosity measurements at LHC with ATLAS was proposed in [2,3]. The 
detail study of the trigger issues and the background conditions has shown 
that using of — > + /i"*"//" with high-pr muons looks favorable. 

Due to the highp-r threshold of LHC experiments (> 1 GeV), the cross section 
of pp — > PP + /x^/x" is small (1 ^ 15 pb) therefore this process can not be used 
for the luminosity monitoring but provides for the absolute calibration of any 
stable high-rate monitor with the statistical accuracy of 1 2 % after of a 
few months operation. 

Below the proposal for the ATLAS offline luminosity determination using the 
two-photon production of //''"//"-pairs [4,7,5] is reported and the preliminary 
signal and background cross section estimates for the LHCb detector are pre- 
sented. 

It should be mentioned that the application of the forward e'^'e^-pairs is un- 
der study in ATLAS [6,7,5]. The possibility to employ pp — > pp -|- /*"*"/*" for 
luminosity measurements was manifested in [8] without a detail consideration 
of theoretical uncertainties and detector related questions. 



2 On possibility of the cross section calculation 

A review of the two-photon leptoproduction can be found elsewhere [9]. Main 
features of this process can be illustrated in the Equivalent Photon Approxi- 
mation (E.P.A.) allowing one to express its cross section via the cross section 
of the pair production by virtual photons with the equivalent photon spectra 

da — (7^*7* -^LL drii dn2 (1) 



For colliding leptons one has [9] 



where q{u!,q) is the four-momentum of the virtual photon and 7 = E/m is 
the Lorentz factor of the coUiding particle (we assume u <^ E). 

The cutoff of the spectrum (2) occurs at Iq^l due to the q^ dependence 

of a^*^*^LL is the invariant mass of the pair produced). The characteristic 
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value of the the total transverse momentum of the pair produced pr ~ {oJi + 
1^2) /t ~ 10 MeV in the luminosity measurements context. 

When pp collide, the proton can dissociate during the photon emission there- 
fore the elastic and the inelastic two-photon processes should be distinguished. 

In the elastic case Fig. (l,a) only the small modification of the equivalent 
photon spectrum is required: 



dUel — dUQED x 



1 — q'^/Arrip 

{Ge, Gm are the electromagnetic form factors). 



(3) 
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(a) (b) (c) 

Fig. 1. Two-photon production of lepton pairs in 'pp collisions: (a) - elastic, (b,c) - 
inelastic 

For the inelastic processes Fig. (l,b), (l,c) one has for each inelastic vertex 
[9,1] 



dn 



inel 



duQED X dM^ ^ duQED X — ^ dM" (4) 
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where M is the invariant mass of the hadronic system, W2 is the inelastic 
scattering structure function and cjj?, are the ^p cross sections known 
from the photo- and electro-production experiments. Unlike (2) and (3), the 
expression (4) is not singular for q'^ — > thus the characteristic pair p^ for the 
inelastic production is not small (~ 250 MeV). 

Diagrams of Fig. 1 reflect the strong interaction inside the single proton, be- 
sides, the strong interaction between colliding protons, so called rescattering, 
should be taken into account (Fig. 2). 

In the ATLAS luminosity measurement context, the rescattering has been 
considered in [10]. It was shown that the elastic strong interaction Fig. (2,b) 
just modifies the phase of the matrix element and does not change the cross 
section. The inelastic strong interaction Fig. (2,c) reduces the yield of pp 
PP + fi~^fi~ on behalf of pp ^ X + In the experimental setup when only 

the muon pair is detected, the total dimuon cross section remains unchanged. 
To reduce the dimuon background from Drell-Yan process and the hadron 
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(a) (b) (c) 

Fig. 2. Modification of the elastic two-photon process by rescattering (final state 
rescattering is not shown). 

decays, the absence of extra particles in the dimuon vertex is necessary. In 
this case 



pp— el I inel /-i \ 

77 ^inel ' ^^77 v-*- ^reseat) ' 



reseat 
77 



(5) 



where tineu ^reseat are the probabilities that the event passes the vertex cut; 
for LHCb and ATLAS eineh (^reseat > 0.8. 

The first two terms can be estimated using the E.P.A. formulae 



da. 



el 

77 



(6) 



cicT^^f' = {dUel,! dninel,2 + dUinel,! dnel,2 + dUinel,! dninel,2) ' O. 



77-»/i+/i- 



(7) 



(7^"^' ~ (7^5y for the pair mass W > \ GeV in absence of cuts on the total pair 
Vt- 

The cr^^''"^"* can not be expressed in terms of the equivalent photon spectra. 
Numerically, a;;*'^"* ~ 0.1% of a^^^ for pT(pair) < 30 MeV [10]. 

The M.C. calculations performed by us for the ATLAS conditions has demon- 
strated that al^^'^"-^ ~ 0.1 (7^!!jf' and its dependence of on the pair is similar 
to that of (J^!!|f . So, it causes no additional problems for the luminosity 
measurements. Below we include the rescattering effect in the inelastic cross 
section. 

It should be noted that the rescattering for the inelastic diagrams (l,b), (l,c) 
is not considered in [10] and is ignored in our cross section calcTilations. This 
and the bad knowledges of the soft inelastic scattering structure functions 
make the inelastic cross section estimates rather uncertain. 

For the forward e"'"e~-pair production the smallness of the pair mass W ensures 
the smallness of and makes the inelastic contribution negligible as was stated 
in [1]. The situation is not so simple concerning the muon pair production 
with pr(A*) > 6 GeV and W > 12 GeV. In this case the inelastic background 
is not negligible even with severe experimental cuts (see Fig. 3). However, the 
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distributions of the pair Pt and related variables are very different allowing 
one to rcdnce the background to an acceptable level and then subtract it as 
described below. 

The uncertainty in the residual inelastic contribution is about 30%, ~ 3% 
of the visible signal, does not affect the luminosity measurement due to the 
background subtraction. The accuracy of the clastic signal calculation < 0.5% 
is limited by the accuracy of the proton form factors only. 

The formulae (6), (7) are presented above just for the illustration. The actual 
calculations were performed using the event generator [11] implementing the 
exact lowest order matrix element. The muon bremsstrahlung, affecting the 
essential parameter distributions, was simulated in the soft photon approxima- 
tion; other types of radiative corrections were not applied. It seems sufficient 
for the current stage of the work. 



3 Procedure for signal extraction 



The following criteria are proposed for the ATLAS event selection: 

(1) Two muon tracks of opposite charges (measured in both the Inner De- 
tector and the muon spectrometer and triggered by the Muon systems) 
with Pt > 6 GeV and I77I < 2.2; 

(2) Muon pair invariant mass < 60 GeV (against Z°); 

(3) Pt of the muons are equal within 2.5(7 of the measurement uncertainty 
{a{pT)/pT ^ 1.5% for Pt < 20 GeV); 

(4) AcoUinearity angle © > 1° so the muons are not exactly back-to-back 
(against the cosmic ray background); 

da/dp^ (pb/MeV) 
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0.02 



0.01 
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elastic yy with 
detector smearing 
inelastic yy 
Drell-Yan with the 
vertex multiplicity cut 
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10, V~ pair Pt, (MeV) 

Fig. 3. The differential dimuon cross section da/dpT for pp ^ pp + ji^ as a 
function of the pair transverse momentum pT (ATLAS, W > 12 GeV, \rf\ < 2.3). 
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(5) Probability of > 1% for the muon vertex fit; 

(6) No other charged tracks with pt > 0.5 GeV, \r]\ < 2.5. (j{z.^,) < 1 mm and 
a good from the muon vertex {a{zy) is the estimated error in z for the 
nearest point to the beam). 

The condition (6) strongly reduces the background from Drell Yan process 
and hadron decays but makes the detection efficiency dependent on the event 
pile-up probability. For the longitudinal LHC bunch size of 7.5 cm and the 
luminosity L — 10^^ cm~^ s~^, the loss of two-photon events is about 3%. 

The expected acoplanarity distribution for selected events in shown in Fig. 4. 
(the acoplanarity angle (p is defined as the angle between the muon production 
planes). The simulation was performed using the particle level Monte Carlo 
code [4] with the detector properties parameterized according to the ATLAS 
specification. The signal and background cross section estimates for the refer- 
ence region |0| < 5 mrad containing about 50% of the signal are presented in 
Table 1. 

Ao/A(|) (pb/0.2 mrad) 
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Fig. 4. The acoplanarity distributions for the signal and dominant background pro- 
cesses with all cuts applied. 

The acoplanarity distribution has the sharp signal peak at = and is almost 
uniform for background processes in the range of |0| < 50 mrad. Unlike the 
distribution of the pair transverse momentum, it is not much affected by the 
detector resolution. Using this, the final signal extraction can be performed by 
fitting the (f) distribution in the range of ±20-=-40 mrad with the signal shape 
obtained by Monte Carlo and the symmetric parabohc background A — BcfP'. 

With the proper modification of the event selection criteria and the back- 
ground subtraction (f) interval, the procedure describe above can be employed 
in CMS and LHCb experiments. 
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Table 1 

ATLAS cross section estimates for the kinematic cuts (1-4) and all cuts (the ex- 
pected detection efficiency > 0.65 is not taken into account). 



process 


<ykin, pb 


<^alh pb 


signal 


1.33 


1.30" 


77 (inelastic) 


0.13 


0.12'' 


Drell-Yan process 


4.0 


0.04 


Heavy quark decays 


10. 


0.01 


tt/K decays 


1.8 


<0.001 


background 


15.9 


0.17 (13%) 



a) for L = 10^^ cm ^ s ^ 

b) including rescattering effect 



4 On accuracy of ATLAS luminosity determination 

The expected statistical accuracy of the luminosity determination ^ 1.5% for 
the integrated luminosity of 10 fb~^ (about 4 months of operation at L = 10^^ 
cm~^s~^) and the detection (trigger and reconstruction) efficiency of 0.65. 

The systematic errors due to the uncertainty in the proton form factors and 
the background subtraction procedure are expected to be < 0.5%. The pile- 
up correction can be done using the experimental data. The systematic error 
should not exceed 1% at L = 10^^ and is neghgible at L = 10^^ cm~^ s~^. 

The uncertainty of the detection efficiency seems dominating source of the 
systematic error. This efficiency can be determined experimentally using the 
pp ^ pp + fj.^fi~ event sample with the accuracy better than 1% (it can be 
keep clean enough with some selection criteria relaxed) . For the low luminosity 
operation {L = 10^^ cm~^s~^) the ffist level trigger on the single muon with 
Pt > Q GeV is foreseen [12], and the trigger efficiency can be determined using 
the second muon (a minor modification of the second level trigger software is 
required). For the operation at L = 10^^ cm~^ the single muon trigger with 
Pt > 6 GeV will have too high rate therefore two such muons will be required 
in the ffist level trigger. The trigger efficiency for this case is supposed to 
be known from the low luminosity run. The other option is triggering on the 
common dimuon track in (p, (^)-projection as shown schematically in Fig. 5. 
This allows additionally to reduce the pr threshold and substantially increase 
the signal rate. This needs some extending of the ATLAS first level trigger. 

Without the experimentally measurement efficiency one has to rely on the 
Monte Carlo calculations involving tracking of low momentum muons in the 
detector material and in the inhomogcueous magnetic field of the air core 
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Fig. 5. The common fi'^/j. track in (p, (/))-projection (the multiple scattering is not 
taken into account). 

toroid. In this case the accuracy of the luminosity measurement with pp — > 
pp + ii'^ii" would hardly be better than 5%. 



5 Cross section estimates for LHCb 



To study the potential of the method suggested for the LHCb experiment, the 
cross section estimates has been obtained using the particle level Monte Carlo 
code. The particles were tracked through some model of the LHCb detector 
without interaction with the material. The resolution in and 9 angles was 
fixed by the number and the width of the strips hitted in the vertex detector. 
The total momentum resolution of 0.3 p % was assumed [13]. 

The following event selection criteria were applied: 

(1) Two muon tracks of opposite charges with pt > 1 GeV, E > S GeV are 

detected; 

(2) Muon pair invariant mass < 30 GeV; 

(3) Difference in the measured values of muon pt is small: 

IPt-PtIUPt+Pt) <0.13. 

(4) No other charged tracks are observed in the event. 

The preliminary results obtained for the reference region |0| < 25 mrad are 
presented in Table 2. The observed multiplicity for all background sources 
except the two-photon one are 30 60. No events passing the condition (4) 
were generated, the suppression factors were roughly estimated using the fits 
of the multiplicity distributions for the relaxed cut (3). 

Due to lower pr threshold, the estimated LHCb signal cross section is more 
than 10 times larger then that for ATLAS. The statistical accuracy of 1% 
seems reachable at the designed luminosity of 2.10^^ cm~^s~^. 

By the same reason the variation of the background level in the background 
subtraction interval |0| < 100 mrad is bigger than for ATLAS. Some opti- 
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Table 2 

LHCb cross section estimates for the kinematic cuts (1-3) and all cuts (the detection 
efficiency is not taken into account). 



process 


CTkin, pb 


<^all, pb 


signal 


14.6 


Ril4.6 


77 (inelastic) 


0.6 


0.5 


Drell-Yan process 


93 


<0.02 


Heavy quarks decays 


690 


<0.02 


TT, K decays 


~50000 


<0.3 


background 


93.6 


0.5 ^ 0.84 (<6%) 



mization of cuts might be necessary to keep the systematic error smalL 
6 Conclusion 

The method described allows for the ATLAS offline luminosity determination 
using the two-photon process pp ^ pp + ii'^ ii~ with the accuracy < 2% in the 
luminosity range of 10^^-^10^^ cm~^s~^. The preliminary cross section esti- 
mates done for LHCb promise the same level of the luminosity measurement 
accuracy at L = 2 • 10^^ cm~^ s~^. 
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